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Трансформации региональньх аномалий поля силь тяжести «(х) во внешнем прост- 


ранстве у с адекватной точностью целесообразно осуществлять, решая нелинейную гранич- 


ную задачу Алексидзе |! | для уравнения Лапласа, в граничньх условиях которой стоят сами 
значения сильг тяжести, являющиеся значениями модуля градиента потенциала сильт тяжести 


(МГПСТ) «(х)е ягад (ж) . Обобщенно постановка задачи Алексидзе такова: найти потенци- 
ал У) х є у", удовлетворяющий внутри замкнутой области у - у бу уравнению Лапла- 
са ДУ (х) -0, хє у", авточках границь ду Ляпунова области у" и наоо - условиям: 


з(ауоо) 
2; жо 2-8 (), х, єду 
ке Ху , 


та УУ(х)-» 0 


хо 


() 


- з 
где є(х) - задана, у - ограниченная область пространства К" с тяготеющими массами, у" 
- ее неограниченноге дополнение без тяготеющих масе, ду - поверхность измерений - грани- 
ца областей у и у/. 


Бе решение, в силу избранной однородной слоистой модели средь, определяет по- 
тенциал простого слоя с искомой плотностью, распространенной на контактной поверхности 
типа Ляпунова. Зта плотность определяєтся как решение уравнения силь тяжести, аналити- 
ческое вьтражение которого являєтся следствиєм свойств МГПСТ огад УУ(х) , а конкретньмй 





вид зависит от избранной модели Земли. В плоском случає зта задача решена итерационно. 
Точность ее решения существенно зависит от обусловленности задачи, параметров метода 
(числа и расположения фундаментальньїжх решений), приближений направления искомого 
градиента, направления внешней нормали, вдоль которой вьтчисляют производньєе в гранич- 
ньх условиях (учет векторной природь поля). Однако, ее решение на простьх моделях сре- 
дьт при определенньгх ограничениях на решение - успешно. 

Задача (1) ориентирована на аналитическое продолжениєе поля вне источников анома- 
льньїх масс при заданной с некой точностью модели средьї и погрешностях измерений. На 
параметрьт средьт наложено ограничение в виде кусочной непрерьвности и дифференцируе- 
мости по Фреше ради изображениє приращений решений задачи линейной комбинацией 
приращений параметров модели (вьщделение главной линейной части искомого решения). 

Основное априорное предположениєе о поведений искомого решения обратной задачи 
(линейная независимость оператора решения и его производньх) обеспечиваєтся одним из 
методов |2). Указанньх ограничений достаточно для численной сходимости к точному 
решению. Решение задачи неоднозначно из-за конечного вертикального разрешения, неадек- 
ватного вьтбор модели задачи и/или начальньх приближений решения, что свойственно при- 
роде обратньтх задач с производньми в граничньтх условиях. Точность дискретизации задачи 
зависит и от способа вьтчислений производньхх. 

Нюансь6 численного моделирования таковьт: осцилляция интегральньх ядер операто- 
ров обратной задачи и существенная зависимость их визуализации от избранного способа 
гриддинга. Нивелировать осцилляцию можно, задействовав конструкции типа интеграла 
Шварца для полосьт |3|, а нюансьт графической визуализации нуждаются в зксперименталь- 


ном исследований и вьборе адекватного инструментария. 
При всей сложности численной реализации задачи Алексидзе, очертим круг ее приме- 
нений: уточнение фигурьт Земли и аналитические трансформации гравимагнитньтх полей. 
Остановимся детальнее на восстановлений магнитного потенциала с помощью задачи 
Алексидзе. Если магнитньй потенциал изобразить как сумму У(Х) - (х)ч т(х) нормального 


С(х) для данной зпохи и аномального т(х) магнитньях потенциалов, то восстановить маг- 


нитньй потенциал УС) по заданньм значениям модуля его градиента Х(х)- гад УС) 5 


хе 
хе бу упрощенно можно по алгоритму, в котором последовательно вьютмчисляєм: 
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2. направляющие косинусьт сопло Х, ) з уч )/7, (); 

3. следующее приближение модуля градиента магнитного потенциала 


3 
УДИ () За (2- У Сол )- Хо ХДУ; 
Кеї 
4. дежурное приближение аномального магнитного потенциала: 
Т..(х) у 7.2 (х)соз(п. т) у(х)сов(у, т), хе бу ; (2) 
5. следующее приближение магнитного потенциала: у, )- С (х ун т.(х ). 


В качестве нормального потенциала целесообразно взять аномальное поле однородно 
наманиченного шара: аномальньтй потенциал т(х) отразит поле реальньаох обьектов, неучтен- 


ньх в зтой модели. Кроме данньх 7(х), хєду (по-сути, ДТ, ) известнь уравнение физичес- 


кой поверхности ду Земли, вектор единичной нормали т(х)- со8(х,,т,), К - 3 к ду почти 


- - - я сі 
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0 (х) є Сх, проходящей через точку х, и со5(у, т) - У со5(у, х, |со5(х,, т)), З - (втад У(х). 
ке Ух 
Если поверхность Земли бу - множество Ляпунова, магнитньюшй потенциал У(х), хєу" 


данной зпохи восстанавливаєтся однозначно по напряженности 7 (х), хєду магнитного по- 
ля зтой зпохи. 


Так, на і-м шаге определим со5( п, ,т) - У, )сок(ху т) и граничное условиє 


Р РН у(х)соз(у, т) (3) 
для определения і-го приближения плотности 5) потенциала простого слоя т) из ли- 
нейного интегрального ше а 
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Где и-х- 8. Далее вьічисляєм приближениє 


и і-1-е приближение направляющих косинусов: соб(п,, Х, ) щ 


Последовательность т (х)) из (4) сходится в себе, как следствие теоремьт: если квад- 
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рат є) отношения модулей градиентов аномального и нормального потенциалов пренеб- 
режимо мал по сравнению с «(х), то последовательность решений т (х)) граничньгх задач 


дТ,(х) ре Е.О, хеду (5) 
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(х)-0, хє у", Ши т.(х)-» 0, 
хо» 
сходится к аномальному потенциалу Т(х) хє у при Т(х) са 00.) п 
Впрочем, последнее условие необязательно, ибо граничное условиєе б - Ех), 


хєду проще, чем условиє (3). Для сходимости последовательности т) достаточно, 


Ї 
чтобь 7 (х) « уУ(х) /сов(у, т), хєду. Зта схема восстановления магнитного потенциала т (х)) 


зкономичнее по обьему вьгчислений, чем схема |4Ї. 

Качество численного решения (среднеквадратичное отклонение параметров модели 
задачи (5) и входньгїх данньх от начального приближения) на текущей итерации контролиру- 
ется по критерию невязки при дополнительном условиий минимума разрешающей способнос- 
ти. Оно зависит от мерьт обусловленности задачи и вьгбора начальньх приближений реше- 
ния. Все решения получень для модели шара и нуждаются в проверке на полевьтх данньх. 
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